A Semiclassical Framework for Mixed Quantum Classical Dynamics.

Semiclassical (SC) approximations for quantum dynamic simulations in complex chemical systems range from rigorously accurate methods that are computationally expensive to methods that exhibit near-classical scaling with system size but are limited in their ability to describe quantum effects. In practical studies of high-dimensional reactions, neither extreme is the best choice: frequently a high-level quantum mechanical description is only required for a handful of modes, while the majority of environment modes that do not play a key role in the reactive event of interest are well served with a lower level of theory. In this feature, we introduce modified Filinov filtration as a powerful tool to construct mixed quantum-classical SC theories where different subsystems can be quantized to different extents without introducing ad hoc intersubsystem interaction terms. We demonstrate that these Filinov-based SC methods can systematically tune between quantum and classical limit SC behavior, offering a practical way forward to accurate and computationally efficient simulations of high-dimensional quantum processes.

Real-time dynamics of strongly correlated fermions using auxiliary field quantum Monte Carlo.

Spurred by recent technological advances, there is a growing demand for computational methods that can accurately predict the dynamics of correlated electrons. Such methods can provide much-needed theoretical insights into the electron dynamics probed via time-resolved spectroscopy experiments and observed in non-equilibrium ultracold atom experiments. In this article, we develop and benchmark a numerically exact Auxiliary Field Quantum Monte Carlo (AFQMC) method for modeling the dynamics of correlated electrons in real time. AFQMC has become a powerful method for predicting the ground state and finite temperature properties of strongly correlated systems mostly by employing constraints to control the sign problem. Our initial goal in this work is to determine how well AFQMC generalizes to real-time electron dynamics problems without constraints. By modeling the repulsive Hubbard model on different lattices and with differing initial electronic configurations, we show that real-time AFQMC is capable of accurately capturing long-lived electronic coherences beyond the reach of mean field techniques. While the times to which we can meaningfully model decrease with increasing correlation strength and system size as a result of the exponential growth of the dynamical phase problem, we show that our technique can model the short-time behavior of strongly correlated systems to very high accuracy. Crucially, we find that importance sampling, combined with a novel adaptive active space sampling technique, can substantially lengthen the times to which we can simulate. These results establish real-time AFQMC as a viable technique for modeling the dynamics of correlated electron systems and serve as a basis for future sampling advances that will further mitigate the dynamical phase problem.

Semiclassical dynamics in the mixed quantum-classical limit.

The semiclassical double Herman-Kluk initial value representation is an accurate approach to computing quantum real time correlation functions, but its applications are limited by the need to evaluate an oscillatory integral. In previous work, we have shown that this "sign problem" can be mitigated using the modified Filinov filtration technique to control the extent to which individual modes of the system contribute to the overall phase of the integrand. Here, we follow this idea to a logical conclusion: we analytically derive a general expression for the mixed quantum-classical limit of the semiclassical correlation function-analytical mixed quantum-classical-initial value representation (AMQC-IVR), where the phase contributions from the "classical" modes of the system are filtered while the "quantum" modes are treated in the full semiclassical limit. We numerically demonstrate the accuracy and efficiency of the AMQC-IVR formulation in calculations of quantum correlation functions and reaction rates using three model systems with varied coupling strengths between the classical and quantum subsystems. We also introduce a separable prefactor approximation that further reduces computational cost but is only accurate in the limit of weak coupling between the quantum and classical subsystems.

Nonadiabatic semiclassical dynamics in the mixed quantum-classical initial value representation.

We extend the Mixed Quantum-Classical Initial Value Representation (MQC-IVR), a semiclassical method for computing real-time correlation functions, to electronically nonadiabatic systems using the Meyer-Miller-Stock-Thoss (MMST) Hamiltonian in order to treat electronic and nuclear degrees of freedom (dofs) within a consistent dynamic framework. We introduce an efficient symplectic integration scheme, the MInt algorithm, for numerical time evolution of the phase space variables and monodromy matrix under the non-separable MMST Hamiltonian. We then calculate the probability of transmission through a curve crossing in model two-level systems and show that MQC-IVR reproduces quantum-limit semiclassical results in good agreement with exact quantum methods in one limit, and in the other limit yields results that are in keeping with classical limit semiclassical methods like linearized IVR. Finally, exploiting the ability of the MQC-IVR to quantize different dofs to different extents, we present a detailed study of the extents to which quantizing the nuclear and electronic dofs improves numerical convergence properties without significant loss of accuracy.

Validating and implementing modified Filinov phase filtration in semiclassical dynamics.

The Mixed Quantum-Classical Initial Value Representation (MQC-IVR) is a recently introduced approximate semiclassical (SC) method for the calculation of real-time quantum correlation functions. MQC-IVR employs a modified Filinov filtration (MFF) scheme to control the overall phase of the SC integrand, extending the applicability of SC methods to complex systems while retaining their ability to accurately describe quantum coherence effects. Here, we address questions regarding the effectiveness of the MFF scheme in combination with SC dynamics. Previous work showed that this filtering scheme is of limited utility in the context of semiclassical wavepacket propagation, but we find that the MFF is extraordinarily powerful in the context of correlation functions. By examining trajectory phase and amplitude contributions to the real-time SC correlation function in a model system, we clearly demonstrate that the MFF serves to reduce noise by damping amplitude only in regions of highly oscillatory phase leading to a reduction in computational effort while retaining accuracy. Further, we introduce a novel and efficient MQC-IVR formulation that allows for linear scaling in computational cost with the total simulation length, a significant improvement over the more-than quadratic scaling exhibited by the original method.

Thermodynamics of peptide dimer formation.

The Replica Exchange Statistical Temperature Molecular Dynamics algorithm is used to study the equilibrium properties of a peptide monomer and dimer and the thermodynamics of peptide dimer formation. The simulation data are analyzed by the Statistical Temperature Weighted Histogram Analysis Method. Each 10-residue peptide is represented by a coarse-grained model with hydrophobic side chains and has an α-helix as its minimum energy configuration. It is shown that the configurational behavior of the dimer can be divided into four regions as the temperature increases: two folded peptides; one folded and one unfolded peptide; two unfolded peptides; and two spatially separated peptides. Two important phenomena are discussed: in the dimer, one peptide unfolds at a lower temperature than the isolated monomer and the other peptide unfolds at a higher temperature than the isolated monomer. In addition, in the temperature region where one peptide is folded and the other unfolded, the unfolded peptide adopts an extended structure that minimizes the overall surface area of the aggregate. It is suggested that combination of destabilization due to aggregation and the resulting extended configuration of the destabilized peptide could have implications for nucleating ß-sheet structures and the ultimate formation of fibrils.